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(57) ABSTRACT 

Mapping coherent/incoherent acoustic sources as determined 
from a phased microphone array. A linear configuration of 
equations and unknowns are formed by accounting for a 
reciprocal influence of one or more cross-beamforming char- 
acteristics thereof at varying grid locations among the plural- 
ity of grid locations. An equation derived from the linear 
configuration of equations and unknowns can then be itera- 
tively determined. The equation can be attained by the solu- 
tion requirement of a constraint equivalent to the physical 
assumption that the coherent sources have only in phase 
coherence. The size of the problem may then be reduced using 
zoning methods. An optimized noise source distribution is 
then generated over an identified aeroacoustic source region 
associated with a phased microphone array (microphones 
arranged in an optimized grid pattern including a plurality of 
grid locations) in order to compile an output presentation 
thereof, thereby removing beamforming characteristics from 
the resulting output presentation. 
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DECONVOLUTION METHODS AND 
SYSTEMS FOR THE MAPPING OF 
ACOUSTIC SOURCES FROM PHASED 
MICROPHONE ARRAYS 

5 

RELATED APPLICATIONS 

This application is a continuation-in-part of the pending 
application Ser. No. 11/126,518, filed May 10, 2005 now U.S. 
Pat. No. 7,783,060 and claims priority to provisional patent 
application Ser. No. 60/914,451 filed on Apr. 27, 2007. 

ORIGIN OF THE INVENTION 

This invention was made by employees of the United States 
Government and may be manufactured and used by or for the 
Government of the United States of America for governmen- 
tal purposes without the payment of any royalties thereon or 
therefor. 20 

TECHNICAL FIELD 

Embodiments are generally related to phased microphone 
arrays. Embodiments are also related to devices and compo- 25 
nents utilized in wind tunnel and aeroacoustic testing. 
Embodiments additionally relate to aeroacoustic tools uti- 
lized for airframe noise calculations. Embodiments also 
relate to any vehicle or equipment, either stationary or in 
motion, where noise location and intensity are desired to be 30 
determined. 

BACKGROUND OF THE INVENTION 

The specification of pending patent application Ser. No. 35 
11/126,518, filed May 10, 2005, is hereby incorporated by 
reference in its entirety for its teaching (herein referred to as 
“the referenced Ser. No. 1 1/126,518”). 

Wind tunnel tests can be conducted utilizing phased micro- 
phone arrays. A phased microphone array is typically config- 40 
ured as a group of microphones arranged in an optimized 
pattern. The signals from each microphone can be sampled 
and then processed in the frequency domain. The relative 
phase differences seen at each microphone determines where 
noise sources are located. The amplification capability of the 45 
array allows detection of noise sources well below the back- 
ground noise level. This makes microphone arrays particu- 
larly useful for wind tunnel evaluations of airframe noise 
since, in most cases, the noise produced by wings, flaps, struts 
and landing gear models will be lower than that of the wind 50 
tunnel environment. 

The use of phased arrays of microphones in the study of 
aeroacoustic sources has increased significantly in recent 
years, particularly since the mid 1990’s. The popularity of 
phased arrays is due in large part to the apparent clarity of 55 
array -processed results, which can reveal noise source distri- 
butions associated with, for example, wind tunnel models, 
and full-scale aircraft. Properly utilized, such arrays are pow- 
erful tools that can extract noise source radiation information 
in circumstances where other measurement techniques may 60 
fail. Presentations of array measurements of aeroacoustic 
noise sources, however, can lend themselves to a great deal of 
uncertainty during interpretation. Proper interpretation 
requires knowledge of the principles of phased arrays and 
processing methodology. Even then, because of the complex- 65 
ity, misinterpretations of actual source distributions (and sub- 
sequent misdirection of engineering efforts) are highly likely. 


2 

Prior to the mid 1980’s, processing of array microphone 
signals as a result of aeroacoustic studies involved time delay 
shifting of signals and summing in order to strengthen con- 
tributions from, and thus “focus” on, chosen locations over 
surfaces or positions in the flow field. Over the years, with 
great advances in computers, this basic “delay and sum” 
processing approach has been replaced by “classical beam- 
forming” approaches involving spectral processing to form 
cross spectral matrices (CSM) and phase shifting using 
increasingly large array element numbers. Such advances 
have greatly increased productivity and processing flexibility, 
but have not changed at all the interpretation complexity of 
the processed array results. 

Some aeroacoustic testing has involved the goal of forming 
a quantitative definition of different airframe noise sources 
spectra and directivity. Such a goal has been achieved with 
arrays in a rather straight-forward manner for the localized 
intense source of flap edge noise. For precise source localiza- 
tion, however, Coherent Output Power (COP) methods can be 
utilized by incorporating unsteady surface pressure measure- 
ments along with the array. Quantitative measurements for 
distributed sources of slat noise have been achieved utilizing 
an array and specially tailored weighting functions that 
matched array beam patterns with knowledge of the line 
source type distribution for slat noise. Similar measurements 
for distributed trailing edge noise and leading edge noise 
(e.g., due in this case to grit boundary layer tripping) have bee 
performed along with special COP methodologies involving 
microphone groups. 

The deconvolution methodology described in the refer- 
enced Ser. No. 1 1/126,518 gives a unique robust deconvolu- 
tion approach designed to determine the “true” noise source 
distribution over an aeroacoustic source region to replace the 
“classical beamformed” distributions. However, that method, 
along with classical beamforming processing, employs sta- 
tistically independent (incoherent) noise source distribution 
assumptions. Thus, it can produce results that are inaccurate 
and distorted in the presence of coherent sources, albeit a 
suitable solution for where non-coherent sources are 
involved. Using an equation form similar to that employed in 
the referenced Ser. No. 1 1/126,518, a solution appropriate to 
identify and quantify coherent as well as an incoherent 
sources is viable and will be herein fully described. 

Example applications for the present invention include 
ideal point and line noise source cases, as well as conforma- 
tion with well documented experimental airframe noise stud- 
ies of wing trailing and leading edge noise, slat noise, and flap 
edge/flap cove noise. 

BRIEF SUMMARY 

The following summary is provided to facilitate an under- 
standing of some of the innovative features unique to the 
embodiments disclosed and is not intended to be a full 
description. A full appreciation of the various aspects of the 
embodiments can be gained by taking the entire specification, 
claims, drawings, and abstract as a whole. 

It is, therefore, one aspect of the present invention to pro- 
vide for a method and system for mapping acoustic sources 
determined from microphone arrays. 

It is another aspect of the present invention to provide for a 
“Deconvolution Approach for the Mapping of Acoustic 
Sources” (DAMAS) when such sources are coherent as well 
as incoherent (DAMAS-C), as determined from phased 
microphone arrays. 
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It is yet a further aspect of the present invention to provide 
for improved devices and components utilized in wind tunnel 
and aeroacoustic testing. 

It is also an aspect of the present invention to provide for 
aeroacoustic tools utilized for airframe noise calculations. 5 

The aforementioned aspects and other objectives and 
advantages can now be achieved as described herein. A 
method and system for mapping coherent and incoherent 
acoustic sources determined from a phased microphone array, 
comprising a plurality of microphones arranged in an opti- to 
mized grid pattern including a plurality of grid locations 
thereof. Utilizing a method similar to that employed in the 
referenced Ser. No. 11/126,518, a linear configuration of 
equations and unknowns can be formed. The present method 
differs in that the terms of the equation are complex and the 15 
problem size for the same number of grid points is expanded. 
The DAMAS-C problem contains N(N+ 1 )/2 potentially inde- 
pendent equations and unknowns. Certain methods are used 
to reduce the computational requirements of solving such a 
system. One or more equations among the linear configura- 20 
tion of equations and unknowns can then be iteratively deter- 
mined. 

In the referenced Ser. No. 11/126,518, the full-rank was 
attained by the solution requirement of the positivity con- 
straint equivalent to the physical assumption of statically 25 
independent noise sources at each location. In the present 
application a similar restriction assumption is made where the 
coherence solutions should be specifically phase related (for 
applications of DAMAS-C where sources are limited as hav- 
ing only in-phase coherence, the constraint sets the value of 30 
the result from a previous iteration to zero when that result is 
not positive). Due to the significant computational require- 
ments of applying DAMAS-C, a further reduction via zoning 
is employed Zoning is a method whereby evaluation is 
restricted to the possible solutions to anticipated or realizable 35 
conditions of the noise source evaluation region under study. 

A noise source distribution is then generated over identified 
aeroacoustic source regions associated with the phased 
microphone array in order to compile an output presentation 
thereof, in response to iteratively determining at least one 40 
equation among the linear configuration of equations and 
unknowns. 

DESCRIPTION OF THE DRAWINGS 

45 

The accompanying figures, in which like reference numer- 
als refer to identical or functionally- similar elements 
throughout the separate views and which are incorporated in 
and form a part of the specification, further illustrate the 
embodiments and, together with the detailed description, 50 
serve to explain the embodiments disclosed herein. 

FIGS. 1A-D illustrate the output dB level contours over 
scan planes of beamforming and crass-beamforming for a 
single source. 

FIG. 2 illustrates a stack of individual n 0 planes defining a 55 
survey and solution space. 

FIGS. 3A-D illustrate results of DAMAS source strengths 
and cross strengths between grid points at n 0 and n over scan 
planes corresponding to FIGS. 1A-D. 

FIGS. 4A-F illustrate beamforming and corresponding 60 
results for both DAMAS and DAMAS-C based on two inco- 
herent point sources being evaluated. 

FIGS. 5A-F illustrate beamforming and corresponding 
results for both DAMAS and DAMAS-C based on two coher- 
ent point sources being evaluated. 65 

FIGS. 6A-F illustrate beamforming and corresponding 
results for both DAMAS and DAMAS-C based on two inco- 
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herent point sources located closer together than the sources 
used in FIGS. 4A-F being evaluated. 

FIGS. 7A-F illustrate beamforming and corresponding 
results for both DAMAS and DAMAS-C based on two coher- 
ent point sources located closer together than the sources used 
in FIGS. 5A-F being evaluated. 

FIGS. 8A-F illustrate beamforming and corresponding 
results for both DAMAS and DAMAS-C based on two inco- 
herent simulated line sources comprised of several point 
sources in a line. 

FIGS. 9A-F illustrate beamforming and corresponding 
results for both DAMAS and DAMAS-C based on two coher- 
ent simulated line sources comprised of several point sources 
in a line. 

FIG. 10 illustrates the setup used to conduct a flap noise test 
in a Quiet Flow Facility. 

FIGS. 11A-H illustrate beamforming and corresponding 
results for both DAMAS and DAMAS-C based on data gath- 
ered from the flap noise test. 

FIG. 12 illustrates a block diagram of a system adapted for 
mapping coherent and incoherent acoustic sources deter- 
mined from a phased microphone array. 

FIG. 13 illustrates a flow diagram of a method for mapping 
coherent acoustic sources determined from a phased micro- 
phone array. 

DETAILED DESCRIPTION 

The particular values and configurations discussed in these 
non-limiting examples can be varied and are cited merely to 
illustrate at least one embodiment and are not intended to 
limit the scope thereof. Additionally, acronyms, symbols, and 
subscripts utilized herein are summarized below. 

SYMBOLS AND ACRONYMS 

a m shear layer refraction amplitude correction for e m/1 
A c DAMAS-C matrix with A^ , , ark components 
A„ , , reciprocal influence of cross-beamforming char- 
acteristics between grid points 
B array half-power “beamwidth” of 3 dB down from beam 
peak maximum 

c 0 speed of sound in medium in the absence of mean flow 

CSM cross spectral matrix 

Y« 0 « 2 coherence between sources at r^ and n 

DR diagonal removal of G in array processing 

Q n steering vector for array for focus at grid point n 

Q mn component of e„ for microphone m 

f frequency 

Af frequency bandwidth resolution of spectra 
G mm , cross- spectrum between P m and P OT , 

G matrix (CSM) of cross- spectrum elements c mm , 

H height of chosen scan plane 
i iteration number 

m microphone identity number in array 
m' same as m, but independently varied 
m 0 total number of microphones in array 
n grid point number on scanning plane(s) 
n'^iojDo 1 same as n but independently varied 
M wind tunnel test Mach number 
X total number of grid points over scanning plane(s) 
p m Fourier Transform of pressure time history at micro- 
phone m 

QFF Quiet Flow Facility 

Q„ idealized p m for modeled source at n for quiescent 
acoustic medium 
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r c distance r m for m equal to the center of the microphone 
array 

x m c 0 retarded coordinate distance from focus point to 
SADA Small Aperture Directional Array 
STD standard or classical array processing 
T complex conjugate transpose (superscript) 
x m propagation time from grid point to microphone m 
w m frequency dependent shading (or weighting) for m 
W shading matrix of w m terms 
W width of scanning plane 
Ax widthwise spacing of grid points 
X c matrix of (X wow ) terms 

n nono (auto) spectrum of “noise source” at grid point n 0 
with levels defined at array, Q* wo Q wo 
X non cross- spectrum between sources at n c and 
n( 

Ay heightwise spacing of grid points 
Y c matrix of Y non terms 

Y nono beamform power response of array at focus location 
n Q Y n of ref app. 

Y non cross -beamform power response between locations 
n Q and n 

The first step in a DAMAS-C formulation is to cross beam- 
form over the source region. FIGS. 1A-D illustrate graphs 
representing output dB contours over scan planes of Beam- 
forming, and Y non and cross beamforming. The referenced 
Ser. No. 11/126,518 describes, in detail, traditional beam- 
forming methods. For the present analysis, the cross beam- 
form product is indicated in equation (1) below: 


The cross-spectral matrix (CSM) is G, where 


However, it is presently desired to find a more general 
distribution for the CSM than that of a distribution of uncor- 
related sources at different n. Using equation (5) the cross 
spectrum between microphones m and in' for a distribution of 
5 sources over all N grid points is given by: 


>w=ZZ 


10 


( 6 ) 


This reflects the acoustic pressure perceived at microphone 
m due to the sources at n', i^, n' 0 , and n are generally different 
than that perceived at microphone nf for the same sources . As 
15 in the referenced Ser. No. 1 1/126,518 the G mm , terms of the 
CSM are proportional to the corresponding P* m P m ' terms. 


20 


=ZZ’ 


( 7 ) 


where 

25 X n> ,=Q\, o Q n , (8) 

X„/ w represents the mean-square cross -spectral pressure 
per bandwidth, due to coherent portion between the sources at 
n 0 ’ andn’, at the microphone m including some normalization. 
The value given by X w , n , is the primary objective when using 
30 DAMAS-C. It is important to note that if the sources at n 0 ' and 
n' radiate noise in a statistically independent way then 
X„, w ,=0. Specifically, if X w , w ,=0 when n 0 W the result col- 
lapses to that found in the referenced Ser. No. 1 1/126,518. 

For the case of a coherent source the CSM is G modc with 
35 components given by equation (7) . U sing equation ( 1 ) we find 


Gn G12 ... Gi mQ 

i c 22 : 



(Y^n) 




AY 



Gm o i 


Gm o m o 


and m 0 is the total number of microphones in the array. This 45 
is a beamform cross -spectrum of the array between focused 
locations of grid points at n=n 0 and at another n. The equiva- 
lent steering “vectors” to those in the referenced Ser. No. 

1 1/126,518 are indicated by equations (3) and (4) below: 

50 

^®lPlnQ^2/jQ ■ ■ ■ (5) 


and 


Z Z( g Uw g -Kx 


Where the bracketed term is 



( 9 ) 


( 10 ) 


e„=col[e ln e 2 n ■ ■ ■ e m(yi ] (4) 

Unlike the referenced Ser. No. 1 1/126,518, where presen- 
tations were of beamforming and solutions over the scan 
plane of N points, the present invention often presents results 
over individual r^ planes with grid points n=l , 2, 3, N. FIG. 2 
illustrates all N of the n 0 planes. Notice eachn 0 plane contains 
all cross-beamforming responses Y wo „ over n=l , 2, 3, N which 
includes standard beamform response Y non at n 0 . 

The pressure transform of a microphone is related to a 
modeled source at a position n in the source field by the 
equation as described in the referenced Ser. No. 1 1/126,518 
and by the following equation: 

Pm:n=Qne m :n~ l ( 5 ) 


e m 0 n' o ) e m Q n' J 
55 

Noting that we can look at explicit terms of equation (9) by 
inserting actual values for n'^n^ and n, the following is 
found: 

60 ... 

Y e =AJC c (11) 

Notice equation (1 1) is the same form as used in the refer- 
enced Ser. No. 1 1/126,518. However here X c and Y c have N 2 
complex-number solutions rather than N real-number com- 
65 ponents. A c has N 4 complex-number components rather than 
N 2 real-number components. The components of A c are given 
by: 
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rl n 0 n,n' 0 n' 



( 12 ) 


5 

Where [ ] M , , is defined by equation 1 0, and the order in A c 
is defined by: 


Ac = Xn\ 


Ani 'hi. 12 
A 2 .ll Al2.12 


[ Aw.u 


Ann.nn \ 


( 13 ) 


10 


15 


The above equations contain terms that are complex con- 
jugates of each other. Further, for the diagonal terms of equa- 
tion (1 3) (i .e. when n^mm^n') the value of A c for that element 
is 1 . These relationships explain why in the present formula- 2 q 
tion there are potentially N(N+l)/2 independent equations 
and unknowns. Therefore, taking advantage of the complex 
conjugate relationships the problem is reduced in size from 
that indicated by the equations above. 

It is noteworthy that modified beamforming such as shaded 2 s 
standard, diagonal removal (DR), and shaded DR beamform- 
ing, as described in the referenced Ser. No. 11/126,518 may 
be applied in a similar manner. All such special beamforming 
processes leave the relationships described above equally 
valid. 30 

To begin solving the DAMAS-C inverse problem, consider 
the following component of equation (11): 


8 

The iteration path is consistent with a progression through 
a stack of solution maps. FIG. 2 illustrates N planes, where 
within each, the counting sequence starts in the left bottom 
corner at n=l and increases vertically along each column. For 
each iteration the value of X (z) „ n replaces the previous itera- 
tion X (l_1) value. Thus, equation 16 represents the solution 
to the DAMAS-C inverse problem described by equation 1 1 . 

In the referenced Ser. No. 11/126,518 a positivity con- 
straint was used in order to render the solutions sufficiently 
deterministic. That positivity constraint was physically nec- 
essary. A similar type of constraint is necessary in the present 
application. In DAMAS-C the value of X non is a complex 
quantity of the form Re(X wo „)+Im(X wow ). When n=n G , 
X wow =X wowo is real and positive, this is equivalent to the 
DAMAS X w which is the autospectral pressure- squared (posi- 
tive) amplitude of sources. Thus, for each iteration of equa- 
tion (16), Im(X wowo ) are set to zero and Re(X wowo ) are set to 
zero only if the value is negative. This is equivalent to the 
positivity constraint described in the referenced Ser. No. 
11/126,518. 

When n*n 0 . X wow could be in any of four complex quad- 
rants. In terms of a complex coherence definition this is rep- 
resented by the following equation: 

x non =y n ^/x^/x nono (17) 

Where the coherence is described as: 



y* y 

/y n 0 n' s -n 0 n 

Xn 0 n 0 Xnn 


(18) 


( 19 ) 



(14) 35 


This equation rearranged (with the appropriate special 
relationships noted above accounted for) gives: 


X-rioti ~ hi 0 n 


yi A no n.n' 0 n' Xn' 0 n' + 

n' 0 =l n' = 1 
N N 

2 2 A no nji' 0 n' ^n' 0 n' 

n'= 1 n'=n + 1 


( 15 ) 


45 


This equation is used in an iteration algorithm to obtain the 
source distribution strengths X„„ (or X„ owo ) for all n and cross 
strengths X wow for all combinations of n 0 and n based on the 
following equation. 


y(0 _ y 
A n 0 n ~ 1 n c n 


V V A I ,x (i) + 

/ j / j n n 0 rui. 0 n * n ' 0 n ’ T 
n o = 1 "' =1 
N N 

V V A , ,x (i ~ l) 

/ j / I n 0 n.n 0 n ^ n ’ 0 ri 


( 16 ) 


60 


Notice the similarity of this form to that given in the refer- 
enced Ser. No. 1 1/126,518. The fundamental difference that 65 
arises is that the terms here are complex and the problem size 
for the same number of N grid points is increased. 


Here, is the phase between coherent portions of the 
source at point n with respect to n 0 . Physically, y„ n is inter- 
preted as the coherence factor between the sources at n and n G . 
This can be related to noise emission from unsteady aerody- 
namic related regions over radiating sources or reflections. 

An appropriate constraint based on the above analysis 
could be enforced in the iterations. However, here X wow is 
regarded as an independent variable just like X wow and X ww . 

Regarding the generality of X wo „, there is a remaining 
question about the rank of the DAMAS-C inverse problem. 
This gives rise to concerns about the practicality of solving 
X wow with arbitrary phase. Thus, in the present application 
example sources are limited to those having in-phase coher- 
ence (<D W M =0) For this example, after each iteration in the 
n^n G case, both the imaginary and real parts ofX non are set to 
zero if the value is not already positive. However, the phase 
<D Wo „ can be specified as being functionally or defined as a 
constant. 

Finally, in an effort to manage the large matrices involved 
in evaluation ofY c =A c X c in DAMAS-C applications, reduc- 
tion by zoning is used. Zoning is employed to restrict the 
possible solutions to anticipated conditions of the noise 
source evaluation region under study. The evaluation region 
can be composed of a number of grid point zones, each with 
assumed coherence criteria. The criteria can be uniform over 
the zones or functionally dependent on, for example, the 
point-to-point distance and frequency. 

In the present application, the source evaluation region is 
composed of multiple non-congruent Zones A and B contain- 
ing grid points (n)^ and (n) 5 . Zone A is taken as a region of 
coherent sources, while Zone B is composed of completely 
incoherent sources. This means cross terms X (w , o) ^ (M ^ are 
zero. Similarly, X^,^^ when n 0 ^n. This leads to the zero- 
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ing out of corresponding A c matrix columns. Now to follow 
the iteration scheme, the cross Yc terms and corresponding 
Xc terms and the corresponding matrix rows of A c are elimi- 
nated. Even though the Yc cross terms themselves are not 
zero, their elimination reduces the number of equations and 5 
unknowns while still giving, weight in the solutions to Zone B 
through the auto terms of Y c . The zoning method described 
above can give rise to a substantial reduction in the size of the 
problem, thus making what might otherwise be an untenable 
computational method attractive. 10 

FIGS. 1A-D illustrate graphs representing output dB con- 
tours over scan planes of Beamforming, and Y non and cross 
beamforming Y n between grid points at n G and n. FIGS. 
1A-D actually illustrate the simplest case, with no issues of 
coherence or multiple sources. The example is of a scan plane 1 5 
placed 60 inches from the 7.8 inch diameter SADA micro- 
phone array. The frequency used was 20 kHz. The scan plane 
is comprised of a 15x15 grid pattern of points where 
Ax=Ay=1.5 inches. The half-power beamwidth, B auto , is 
approximately 7 inches. The corresponding beamwidth B cross 20 
is approximately 10.5 inches. Criteria given in the referenced 
Ser. No. 11/126,518 for resolution range is similarly appli- 
cable, and was met here. Those criteria are: 

0.05 (or Ay/i?)=i0.2 (20) 

and 

l^W/B (and H/B). (21) 

All 225 grid points are considered in Zone A, where coher- 30 
ence is permitted. The point source is located at n=l 13. Y n „ 
is then solved over the scan plane and plotted in graph 110 
illustrated in FIGS. 1A-D. 

FIG. 2 illustrates a stack of individual planes 200, repre- 
senting the DAMAS -C results for X wow and X wow correspond- 35 
ing to the Y non ^ andY WoW plots of FIGS. 1A-D. These were 
determined using the algorithm from equation (16) using 
i=100 iterations. FIGS. 3A-D are graphs 300 illustrating lev- 
els Y fl0flii which represent the collection of X^ values, when 
n=n o , from the individual n 0 planes. Graph 310 illustrates that 40 
DAMAS -C, with this number of iterations, approaches the 
correct source definition. The value found for X 113113 was 
96.5 dB, only slightly differing from the actual value of 100 
dB. The nearby X MqM levels were about 89 dB (the exact value 
being -oo). Finally, the result of a total of all the grid points 45 
was 100.2 dB with the exact value being lOOdB. Of course, an 
increase in the number of iterations preformed would 
increase accuracy. Here, the relatively low number of itera- 
tions illustrates “energy” smearing seen when fewer itera- 
tions are preformed. 50 

Next the ability of DAMAS-C to separate and quantify 
different sources was tested. FIGS. 4A-F illustrate diagramed 
results for two incoherent sources positioned 9 inches apart 
on a 51x51 inch scan plane with a grid point spacing of 
Ax=Ay=l inch. The top left frame 410 is an illustration of 55 
beamforming. The top right frame 420 is a DAMAS pro- 
cessed result using methods described in the referenced Ser. 
No. 11/12 6, 518. Since the sources are incoherent, application 
of DAMAS correctly yields location and level, 103 dB (100 
dB being the correct result) of the sources. This frame shows 60 
2000 iterations. 

Subsequent frames of FIGS. 4A-F, in particular frames 430 
and 440, show distributions of X non the sum of X„ n over all 
planes, X„ „ (for source location n o =1097 plane), and X„ „ 
(for source location n 0 =l 556 plane). This shows DAMAS-C 65 
correctly identifies the sources. The results indicate some 
smearing due to the limited number of iterations. Frames 410 
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and 420 illustrate the key result; that DAMAS-C correctly 
separates the sources and validates that the two sources have 
no coherence with one another. 

FIGS. 5A-F illustrate graphical results for two sources 
defined as perfectly coherent and in-phase. The configuration 
used to obtain FIGS. 5A-F is identical to that of FIGS. 4A-F. 
The first frame 510 of FIGS. 5A-F shows a geometrically 
distorted result, although the sum of apparent sources is 
nearly correct. However, DAMAS-C is shown to correctly 
separate and quantify the coherent sources. Although the 
levels are slightly lower than the actual results this is 
explained by the same resolution energy smearing phenom- 
ena as described above. 

FIGS. 6A-F and FIGS. 7A-F are graphs illustrating a simi- 
lar set of results to those illustrated in FIGS. 4A-F and FIGS. 
5A-F. In these experiments the sources were placed closer 
together (spaced by 4.5 inches) with Ax=Ay=0.9 inches. For 
the graphs in FIGS. 6A-F, the sources used were incoherent. 
As expected, almost the same degree of success in source 
definition is found using DAMAS and DAMAS-C. However, 
graphs in FIGS. 7A-F illustrate the difficulty that can be 
encountered when the sources are coherent. In FIGS. 7A-F, 
the totaled YX non frame result appears almost as a line 710 
and the X wow frames show smearing between the sources. 
However, following criteria for resolvability detailed by the 
referenced Ser. No. 11/126,518, the present DAMAS-C 
results appear compatible with these criteria. 

FIGS. 8A-F and FIGS. 9A-F illustrate graphs for a set of 
presentations, similar to the preceding, except that thirteen 
sources are distributed to simulate a 12-inch line source. 
Graphs in FIGS. 8A-F illustrate for the incoherent line 
source, that both DAMAS and DAMAS-C give good spatial 
and level definition. The graphs in FIGS. 9A-F illustrate a 
presentation for a coherent line source. The DAMAS result is 
substantially distorted, although the total levels for both 
DAMAS and DAMAS-C are correct. The results of these and 
the preceding figures validate the correctness and functional- 
ity of the DAMAS-C algorithm. 

DAMAS-C was applied to data collected from an airframe 
noise test in a Quiet Flow Facility. FIG. 10 illustrates the flap 
edge test system 1000 configuration. The SADA array 1010 is 
positioned outside the flow field of the system 1000 at a 
distance of 5 feet from the flow model 1020. In this test the 
flap angle 1030 was set at 29 degrees and M=0.1 1. Sources 
were evaluated along a scanning plane aligned with an airfoil 
main element chordline. 

The DR beamform processing and corresponding DAMAS 
results are shown as graphs in FIGS. 11A-H. Zone A is a 9x24 
point region over the scan plane with grid spacing Ax=Ay=l 
inch. The results illustrated in FIGS. 11A-H illustrate that the 
results for DAMAS and DAMAS-C substantially match in 
source and distribution levels. This suggests the flap edge and 
flap cove noise regions can be regarded as distributions of 
incoherent sources to the extent that is resolvable for this size 
array and processing. 

Referring to FIG. 12, a block diagram 1200 of a system in 
accordance with features of the present invention is illus- 
trated. The system 1200 is adapted for mapping coherent and 
incoherent acoustic sources 1250 determined from a phased 
microphone array and includes a plurality of microphones 
1210-121« arranged in an optimized grid pattern 1230 and 
including a plurality of grid locations thereof. Also included 
in the system 1200 is a computer 1220 connected to the 
plurality of microphones 1210-122«, the computer 1220 
adapted for processing any combination of DAMAS 1222 or 
DAMAS-C 1224 modules including: a linear configuration of 
equations and unknowns formed by accounting for cross- 
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beamforming characteristic thereof at varying grid locations 
among said plurality of grid locations; an equation iteratively 
determined from said linear configuration of equations and 
unknowns based on a DAMAS-C inverse formulation; and an 
optimized noi se source di stribution generated over an identi- 5 
fied aeroacoustic source region 1250 associated with said 
phased microphone array in order to compile an output pre- 
sentation thereof, in response to iteratively determining said 
equation among said linear configuration of equations and 
unknowns. 

Referring to FIG. 13, a flow diagram 1300 of a method for 
mapping coherent acoustic sources determined from a phased 
microphone array that can be followed in accordance with 
carrying out aspects of the present invention is illustrated. As 
shown in Block 131 0, a plurality of microphones arranged in 1 5 
an optimized grip pattern including a plurality of grid location 
thereof and connected to a computer adapted to process 
DAMAS and/or DAMAS-C modules is provided. As shown 
in Block 1320, a linear configuration of equations and 
unknowns are formed using cross-beaming characteristics 20 
thereof at varying grid locations among the plurality of grid 
locations. An equation is then iteratively determined from the 
linear configuration of equations and unknowns based on a 
DAMAS-C inverse formulation, as shown in Block 1330. 
Then, as shown in Block 1340, an optimized noise source 25 
distribution is generated over an identified aeroacoustic 
source region associated with the phased microphone array in 
order to compile an output presentation thereof. Generation 
of the optimized noise source distribution can be in response 
to iteratively determining the equation among the linear con- 30 
figuration of equations and unknowns, thereby removing the 
beamforming characteristic from the output presentation. 

The linear configuration can further comprise a system of 
linear equations including Y^A^X^, wherein said system of 
linear equations relates a spatial field of point locations with 35 
beamformed array-output responses thereof to equivalent 
source distributions at a same location. A variable A among 
the system of linear equations can be utilized to disassociate 
an array thereof from acoustic sources of interest. Solving for 
a variable x among said system of linear equations can further 40 
include the equation Y C =A C X C . The variable X can be allowed 
to be an imaginary number with a real part and imaginary part. 
Iteratively determining the equation among the linear con- 
figuration of equations and unknowns can further include the 
step of attaining the equation utilizing a solution requirement 45 
of a constraint that sets the phase of said variable X. If phase 
is limited to zero then the imaginary part of X is set to zero if 
it is already not positive. Iteratively determining the equation 
among said linear configuration of equations and unknowns 
can also include the step of attaining the equation utilizing a 50 
reduction of the size of the problem by zoning. 

It is important to note that the methodology described 
above with respect to the figures and equations, which is 
referred to generally by the DAMAS or DAMAS-C acronym, 
can be implemented in the context of a module(s). In the 55 
computer programming arts, a module (e.g., a software mod- 
ule) can be implemented as a collection of routines and data 
structures that perform particular tasks or implement a par- 
ticular abstract data type. Modules generally can be com- 
posed of two parts. First, a software module may list the 60 
constants, data types, variable, routines and the like that that 
can be accessed by other modules or routines. Second, a 
software module can be configured as an implementation, 
which can be private (i.e., accessible perhaps only to the 
module), and that contains the source code that actually 65 
implements the routines or subroutines upon which the mod- 
ule is based. 
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Thus, for example, the term “module,” as utilized herein 
generally refers to software modules or implementations 
thereof. The world module can also refer to instruction media 
residing in a computer memory, wherein such instruction 
media are retrievable from the computer memory and pro- 
cessed, for example, via a microprocessor. Such modules can 
be utilized separately or together to form a program product 
that can be implemented through signal-bearing media, 
including transmission media and recordable media. 

Accordingly, a program product for mapping coherent and 
incoherent acoustic sources determined from a phased micro- 
phone array can be provided in accordance with features of 
the present invention. The program product can include a 
plurality of microphones which can be arranged in an opti- 
mized grid pattern including a plurality of grid locations 
thereof, instruction media residing in a computer memory for 
forming a linear configuration of equations and unknowns by 
accounting for a reciprocal influence of a cross-beamforming 
characteristic thereof at varying grid locations among said 
plurality of grid locations, instruction media residing in a 
computer for iteratively determining an equation from said 
linear configuration of equations and unknowns based on a 
DAMAS-C inverse formulation and instruction media resid- 
ing in a computer for generating an optimized noise source 
distribution over an identified aeroacoustic source region 
associated with said phased microphone array in order to 
compile an output presentation thereof, in response to itera- 
tively determining said equation among said linear configu- 
ration of equations and unknowns, thereby removing said 
beamforming characteristic from said output presentation. 
Each of said instruction media residing in a computer can be 
comprised of signal-bearing media. The signal-bearing 
media can also comprise at least one of the following types of 
media: transmission media or recordable media. 

It will be appreciated that variations of the above-disclosed 
and other features and functions, or alternatives thereof, may 
be desirably combined into many other different systems or 
applications. Also that various presently unforeseen or unan- 
ticipated alternatives, modifications, variations or improve- 
ments therein may be subsequently made by those skilled in 
the art which are also intended to be encompassed by the 
following claims. 

What is claimed as new and desired to be secured by 
Letters of Patent of the United States is: 

1. A method for mapping coherent or incoherent acoustic 
sources determined from a phased microphone array, com- 
prising the steps of: 

arranging a plurality of microphones in a grid pattern 
wherein a plurality of grid locations thereof, are defined; 
providing a computer connected to said plurality of micro- 
phones for receiving signals generated by said micro- 
phones in response to sound sensed thereby; 
forming, using said computer, a linear configuration of 
equations and unknowns using cross-beamforming 
characteristics based on said signals generated at vary- 
ing grid locations among said plurality of grid locations; 
iteratively determining, using said computer, an equation 
from said linear configuration of equations and 
unknowns based on a DAMAS-C inverse formulation; 
generating, using said computer, an optimized noise source 
distribution over an identified aeroacoustic source 
region associated with said phased microphone array in 
response to said step of iteratively determining said 
equation among said linear configuration of equations 
and unknowns; and 
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compiling, using said computer, an output presentation of 
said optimized noise source distribution wherein said 
cross-beamforming characteristics are not present 
therein. 

2. The method of claim 1, wherein said linear configuration 
further comprises a system of linear equations comprising 
Y C =A C X C , wherein said system of linear equations relates a 
spatial field of point locations with beamformed array-output 
responses thereof to equivalent source distributions at a same 
location. 

3. The method of claim 2, wherein a variable A among said 
system of linear equations is utilized to disassociate an array 
thereof from acoustic sources of interest. 

4. The method of claim 2, further comprising the step of 
solving for a variable X among said system of linear equa- 
tions comprising Y c =Aji c . 

5. The method of claim 4, further comprising the step of 
allowing said variable X to be an imaginary number with a 
real part and imaginary part. 

6. The method of claim 4, wherein said step of iteratively 
determining said equation among said linear configuration of 
equations and unknowns further comprises the step of attain- 
ing said equation utilizing a solution requirement of a con- 
straint that sets real and imaginary parts of said variable X to 
zero when said parts are not positive after each iteration. 

7. The method of claim 1, wherein said step of iteratively 
determining said equation among said linear configuration of 
equations and unknowns further comprises the step of attain- 
ing said equation utilizing zoning. 

8. A computer system for mapping coherent and incoherent 
acoustic sources determined from a phased microphone array 
formed by a plurality of microphones arranged in a grid 
pattern with a plurality of grid locations thereof being defined 
and said plurality of microphones generating signals in 
response to sound sensed thereby, said 

computer system comprising a computer for forming a 
linear configuration of equations and unknowns by 
accounting for a reciprocal influence of a cross-beam- 
forming characteristic based on said signals at varying 
grid locations among said plurality of grid locations; 

said computer iteratively determining an equation from 
said linear configuration of equations and unknowns 
based on a DAMAS-C inverse formulation; 

said computer generating an optimized noise source distri- 
bution over an identified aeroacoustic source region 
associated with said phased microphone array in 
response to said equation so-iteratively determined 
among said linear configuration of equations and 
unknowns; and 

said computer compiling an output presentation of said 
optimized noise source distribution wherein said cross- 
beamforming characteristic is not present in said output 
presentation. 

9. The computer system of claim 8, wherein said linear 
configuration further comprises a system of linear equations 
comprising Y C =A C X C , wherein said system of linear equa- 
tions relates a spatial field of point locations with beam- 
formed array-output responses thereof to equivalent source 
distributions at a same location. 

10. The computer system of claim 9, wherein a variable A 
among said system of linear equations is utilized to disasso- 
ciate an array thereof from acoustic sources of interest. 
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11. The computer system of claim 9, wherein said com- 
puter solves for a variable^ X among said system of linear 
equations comprising Y C =A Jt c . 

12. The computer system of claim 11, wherein said com- 
5 puter iteratively determining said equation among said linear 

configuration of equations and unknowns further comprises 
attaining said equation utilizing a solution requirement of a 
constraint that sets real and imaginary parts of said variable X 
to zero when said parts are not positive after each iteration, 
to 13. The computer system of claim 8, wherein said com- 
puter iteratively determining said equation among said linear 
configuration of equations and unknowns further comprises 
attaining said equation utilizing zoning. 

14. A system for mapping coherent and incoherent acoustic 
15 sources determinable from phased microphone arrays, com- 
prising: 

a plurality of microphones arranged in a grid pattern 
wherein a phased microphone array is formed with a 
plurality of grid locations thereof being defined, said 
20 plurality of microphones generating signals in response 

to sound sensed thereby; and 
a computer connected to said plurality of microphones, 
said computer: 

processing a linear configuration of equations and 
25 unknowns formed by accounting for cross-beamform- 
ing characteristics based on said signals at varying grid 
locations among said plurality of grid locations, 
iteratively determining an equation from said linear con- 
figuration of equations and unknowns based on a 
30 DAMAS-C inverse formulation, 

generating an optimized noise source distribution over an 
identified aeroacoustic source region associated with 
said phased microphone array in response to said equa- 
tion so-iteratively determined among said linear con- 
35 figuration of equations and unknowns, and 

compiling an output presentation of said optimized noise 
distribution wherein said cross-beamforming character- 
istics are not present in said output presentation. 

15. The system of claim 14, wherein iteratively determin- 
40 ing said equation among said linear configuration of equa- 
tions and unknowns further comprises attaining said equation 
utilizing zoning. 

16. The system of claim 14, wherein said linear configura- 
tion further comprises a system of linear equations relating a 

45 spatial field of point locations with beamformed array-output 
responses thereof to equivalent source distributions at a same 
location. 

17. The system of claim 16, wherein a variable A among 
said system of linear equations is utilized to disassociate an 

50 array thereof from acoustic sources of interest. 

18. The system of claim 16, further comprising solving for 
a variable X among said system of linear equations compris- 
ing y=AX 

19. The system of claim 18, wherein iteratively determin- 
55 ing said equation among said linear configuration of equa- 
tions and unknowns further comprises attaining said equation 
utilizing a solution requirement of a constraint that sets real 
and imaginary parts of said variable X to zero when said parts 
are not positive after each iteration. 

60 



